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INTRODUCTION 


In  a  previous  paper**-,  a  new  finite-difference  approach  was  developed  for 
solving  the  axisymmetric  plane-strain  problems  subjected  to  internal  or  exter¬ 
nal  pressure  beyond  the  elastic  limit.  The  material  was  assumed  to  obey  the 
von  Mises*  yield  criterion,  the  Prandt 1-Reuss  flow  theory  and  the  isotropic 
hardening  rule.  The  ideally-plast ic  material  was  treated  as  a  special  case. 
The  new  formulation  is  simpler  than  other  finite-difference  methods  for 
ideally-plast  ic  materials21  and  strain- hardening  materials  .  The  load  incre¬ 
ments  used  in  all  steps  were  fixed  and  equal.  Accurate  numerical  results  can 
be  obtained  by  reducing  the  grid  sizes*  and  load  increments. 

In  the  present  paper,  an  adaptive  algorithm  to  generate  a  more  accurate 
solution  will  be  developed  for  the  plane-strain  problem  of  an  overstrained 
tube.  The  load  increments  in  all  steps  are  varied  and  determined  automatical¬ 
ly  in  the  program.  One  additional  grid  point  will  become  yielded  in  each  load 
step.  To  reach  100%  overstrain,  the  number  of  steps  are  equal  to  the  number 
of  grids.  For  a  given  percentage  of  overstrain  and  a  desired  solution 
accuracy,  the  stresses  and  strains  can  be  obtained  in  a  much  more  efficient 
manner • 


*Chen,  P.  C.  T.  ,  "A  Finite-Difference  Approach  to  Axisymmetric  Plane-Strain 
Problems  Beyond  the  Elastic  Limit,"  Transactions  Twenty-Fifth  Conference  of 
Army  Mathematicians,  pp.  455-466,  January  1980. 

^Hodge,  P.  G.  and  White,  G.  N. ,  "A  Quantitative  Comparison  of  Flow  and 
Deformation  Theories  of  Plasticity,"  J.  Appl.  Mech. ,  Vol.  17,  1950,  pp. 
180-184. 

^Chu,  S.  C.,  "A  More  Rational  Approach  to  the  Problem  of  an  Elasto-Plast ic 
Thick-Walled  Cylinder,"  J.  of  the  Franklin  Institute,  Vol.  294,  1972,  pp. 
57-65. 
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BASIC  EQUATIONS 

Assuming  small  strain  and  no  body  forces  in  the  axi symmetric  state  of 
plane  strain,  the  radial  and  tangential  stresses,  or  and  o0,  must  satisfy  the 
equilibrium  equation, 

r(3or/3r)  =  Oq  -  or  ;  (1) 

and  the  corresponding  strains,  er  and  e0,  are  given  in  terms  of  the  radial 
displacement,  u,  by 

er  =  3u/3r  ,  e0  =  u/r  .  (2) 

It  follows  that  the  strains  must  satisfy  the  equation  of  compatability 

r ( 3e0/3r)  =  er  -  e0  .  (3) 

The  material  is  assumed  to  be  elastic-plastic,  obeying  Misesf  yield 
criterion,  Prandtl-Reuss  flow  theory  and  isotropic  hardening  law.  The 
complete  stress-strain  relations  can  be  rewritten  in  an  incremental  form^ 


d0i  =  dijdej  for  i,j  =  r,0z  (4) 

and 

dij/20  =  v/(l-2v)  +  6^  -  o^Oj’/S  ,  (5) 

where  E  is  Young’s  modulus,  v  is  Poisson’s  ratio,  is  the  Kronecker  delta, 

2  1  9 

S  =  -  (1  +  -  H’/G)o2  ,  2G  *  E/(l+v)  ,  (6) 

3  3 

°ra  =  (or-H3  0+<Jz)/3  ,  0i»  =  a±  -  om  ,  (7) 

a  =  (l//2)[(ar-a0)2  +  (a0-az)2  +  (aa-ar)2] x/2  >  a0  ,  (8) 


^Yamada,  Y. ,  Yoshimura,  N. ,  and  Sakurni ,  T. ,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  for  the  Solution  of  Elastic-Plastic  Problems  by  the  Finite 
Element  Method,"  Int .  J.  Mech.  Sci.,  Vol.  10,  1968,  pp.  343-354. 
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and  o0  is  the  yield  stress  in  simple  tension  or  compression.  For  a  strain¬ 
hardening  material,  H*  is  the  slope  of  the  effective  stress/plastic  strain 
curve 


a  =  H(JdeP)  . 

For  an  ideally-plast ic  material,  H1  =  0.  When  a  <  oQ  or  da  <  0,  the  state  of 
stress  is  elastic  and  the  third  term  in  equation  (4)  disappears.  Consider  a 
thick-walled  cylinder  of  inner  radius  a  and  external  radius  b.  The  tube  is 
subjected  to  inner  pressure  p  and/or  external  pressure  q.  The  elastic 
solution  for  this  problem  is  well-known 


or 


°0 


■j  A  (q-p)  b  2  p-q(b/a)2 


(9) 


b2/a2-l  'r'  b2/a2-l 

oz  -  v(or+ae)  , 

and  the  pressure  p*  or  q*  required  to  cause  initial  yielding  can  be  determined 
by  using  the  Mises'  yield  criterion. 


FINITE-DIFFERENCE  FORMULATIONS 

For  pressure  beyond  the  elastic  limit,  an  incremental  approach  of  the 
finite-difference  formulation  is  used.  At  the  beginning  of  each  incremental 
loading,  Ap  or  Aq,  the  distribution  of  displacements,  strains  and  stresses  is 
assumed  to  be  known  and  we  want  to  determine  Au,  Aer,  Aeq ,  Aar,  Acjq,  Aaz  at 
all  grid  points.  The  cross  section  of  the  tube  is  divided  into  n  rings  with 

rl  =  a,r2,...,rk=p,...,rn+1=b  (10) 

where  p  is  the  radius  of  the  elastic-plastic  interface. 

Since  the  incremental  stresses  are  related  to  the  incremental  strains  by 
the  incremental  form  (Eq.  (4))  and  Au  =  rAeg,  there  exists  only  two  unknowns 
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at  each  station  that  have  to  be  determined  for  each  increment  of  loading.  The  . 
unknown  variables  in  the  present  formulation  are  (Acq)^,  (Aer)^,  for  i  = 

1 ,  2 , . . . j n  , n+ 1 . 

The  equation  of  equilibrium  (Eq.  (1))  and  the  equation  of  compatibility 
(Eq.  (3))  are  valid  for  both  the  elastic  and  the  plastic  regions  of  a 
thick-walled  tube.  The  finite-difference  forms  of  these  two  equations  at  i  = 
l,...,n  are  given  by 

(ri+1-2ri)(Aar)i  -  (r^-r^ (Aae)1  +  ^(*0,^+! 

=  (ri+1-r1)(ae-or)i  -  r±[ (ar) ±+1  -  (ar)i]  (11) 

for  the  equation  of  equilibrium,  and 

(ri+1-2ri)(Ae6)i  -  (r1+1-rt)  (Aer)  ±  +  r1(Ae0)i+1 

=  (ri+l-ri)(er-e0^i  "  riUe0)i+i  “  (£9)1]  (12) 

for  the  equation  of  compatibility. 

With  the  aid  of  the  incremental  stress-strain  relations  (Eq.  (4)),  Eq.  (11) 
can  be  rewritten  as 

[ (ri+l“2r1) (d12)i  +  (-ri+l+r1)(d22)i] (Aee)i 

+  Kri+l-2ri)(dll)i  +  (-ri+l+ri) (d2l ) il  (Aer)l 
+  r  i  ( d  1 2)  i+1  (  A  Ee  )  i-t-1  +  ri(dll)l+l(Aer^i+l 
=  (ri+1-ri)(ae-ar)i  -  r± [ (or)i+1  -  (ar)i]  .  (13) 

The  boundary  conditions  for  the  problem  are 

Aor(a,t)  =  -Ap  ,  Aor(b,t)  =  -Aq  .  (14) 

Using  the  incremental  relations  (Eq.  (4)),  we  rewrite  Eq.  (11)  as 

(dl2)i(Ae6)i  +  (dn)i(Aer)i  =  -Ap  .  (15) 
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and 


(d12^n+l(Ae6^n+l  +  (dll^n+l(Aer)n+l  "  “A<1  * 

Now  we  can  form  a  system  of  2(n+l)  equations  for  solving  2(n+l)  unknowns, 
(Ae@)i,  (Aer)i,  for  i  -  1 ,2, . . . ,n,n+l .  Equations  (15)  and  (16)  are  taken  as 
the  first  and  last  equations,  respectively,  and  the  other  2n  equations  are  set 
up  at  i  =  1,2,. ..,n  using  Eqs.  (12)  and  (13).  The  final  system  is  an 
unsymmetric  band  matrix  with  the  nonzero  terms  clustered  about  the  main 
diagonal,  two  below  and  one  above. 

INCREMENTAL  LOADING  -  FIXED  VS.  SCALED 

When  the  total  applied  pressure  p  or  q  is  given,  it  is  natural  to  divide 
the  loading  path  into  m  equal  fixed  increments  with 

Ap  =  (p-p*)/m  ,  Aq  -  (q-q*)/m  .  (17) 

These  fixed  increments  need  not  be  equal  for  all  steps  and  any  sequence  of  m 
increments  can  be  supplied  by  the  user.  The  fixed  increments  are  applied 
until  the  total  pressure  or  a  given  percentage  of  overstrain  is  reached.  The 
percentage  of  overstrain  is  defined  as  (p-a)/(b-a)  x  100%.  The  accuracy  of 
the  numerical  results  will  depend  upon  the  values  of  m  and  n  used.  Large 
values  of  m  and  n  will  yield  better  results  at  greater  cost.  For  each  value 
of  n  =  20,  50,  100,...,  we  may  set  m  =  n,  2n,  4n,  8n,...,  to  discuss  the 
convergence.  The  numerical  results  suggest  that  a  sequence  of  decreasing 
load-increments  is  a  better  choice  than  that  of  equal  increments. 
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In  the  following,  a  method  to  generate  a  sequence  of  load-increments  is 
described.  The  method  is  based  on  a  scaled  incremental-loading  approach.1^ 

In  each  step,  a  dummy  load-increment  such  as  Ap  is  applied  and  the  incremental 
results  AOi  for  i  -  r,6,z  at  all  grids  are  determined.  For  all  grid  points  at 
which  o  *  |  |o±|  |  <  o0,  we  compute  the  scaler  a's  by  the  formula 

ot  =-  {r  4-  [T2  +4||Aai||2(a02-||ai||z)]1/2}/||Aai||;2  ,  (18) 

2 

wh.e  re 

r  =  |  | Oil  I2  +  |  1 Aa± | |2  -  I |o±  +  Aa± | |2  ,  (19) 

and  |  |  a  ±  I  |  ,  I  lAaJ  |  ,  I  Ic^  +  AaJ  |2  are  computed  by 

|  | a ^ | | 2  =  -  [(ar-ae)2  +  (a0-az)2  +  (az-ar)2]  .  (20) 

i  2 

Let  X  be  the  minimum  of  the  a's.  Then  X  is  the  load-increment  factor  just 
sufficient  to  yield  one  additional  point.  A  sequence  of  xU)  can  be  deter¬ 
mined  for  all  steps  j  =  1,2 . m  and  the  updated  results  are 

p(j)  =  p(j-l)  +  xU)Ap(j) 

o^(j)  =  +  X^J^Ao^j)  ,  etc.'  (21) 

NUMERICAL  RESULTS  AND  DISCUSSIONS 

The  numerical  results  for  a  thick  tube  with  b/a  =  2,  v  =  .3  subjected  to 
internal  pressure  only  were  obtained.  The  elastic— perfectly-plastic  case  (H 
=  0)  as  well  as  strain-hardening  case  (H'  =  E/9)  were  considered.  Both  fixed 
and  scaled  incremental  loading  approaches  were  used.  The  pressure  p*  required 
to  cause  initial  yielding  is  .4323  aQ. 

4Yamada,  Y.,  Yoshimura,  N. ,  and  Sakurni,  T. ,  "Plastic  Stress-Strain  Matrix  and 
Its  Application  for  the  Solution  of  Elastic-Plastic  Problems  by  the  Finite 
Element  Method,"  Int.  J.  Mech.  Sci.,  Vol.  10,  1968,  pp.  343-354. 
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In  order  to  compare  the  rate  of  convergence  based  on  two  incremental- 
loading  procedures,  we  consider  an  elast ic-perf ectly-plast ic  case.  The 
numerical  results  are  shown  in  Tables  I  and  II.  Based  on  the  equal 
load-increments,  the  effect  of  m  and  n  on  the  internal  pressure  and  bore 
displacement  can  be  seen  in  Table  I.  Larger  values  of  m  and  n  will  yield 
better  results  except  the  bore  displacement  for  100%  overstrain.  If  the 
scaled  incremental-loading  approach  is  used,  the  load-increments  for  all  steps 
are  calculated  and  m  -  n  x  100%  overstrain.  The  numerical  results  for  the 
pressure,  displacement,  axial  stress  at  the  inside  and  the  maximum  hoop  stress 
are  shown  in  Table  II.  A  comparison  of  Tables  I  and  II  indicates  that  the 
scaled  incremental-loading  approach  is  much  more  accurate  and  efficient  than 
the  equal  incremental-loading  approach.  For  example,  to  reach  50%  overstrain 
with  n  =  20,  10  scaled  load-increments  can  give  better  results  than  2044  equal 
load-increment  s. 

Finally,  we  consider  a  strain-hardening  tube  subjected  to  internal 
pressure  only.  The  numerical  results  were  based  on  the  scaled  incremental 
loading  and  the  following  parameters:  b/a  =  2,  v  -  .3,  H1  =  E/9,  n  =  100. 

The  stresses  and  strains  as  functions  of  overstrain  percentage  were  obtained. 
In  Figures  1  and  2,  we  show  the  residual  stresses  resulting  from  50%  and  100% 
overstrain.  The  effect  of  favorable  residual  stresses  of  an  autof rettaged 
tube  is  well-known.  A  simple  and  efficient  approach  to  compute  accurate 
residual  stresses  is  important. 
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TABLE  I.  CONVERGENCE  STUDY  BASED  ON  FIXED  LOAD- INCREMENTS  IN  A  PLANE-STRAIN 
TUBE  (b/a  -  2,  V  -  .3,  H'  =  0) ,  Ap  =  (p-p*)/m 
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TABLE  II.  CONVERGENCE  STUDY  BASED  ON  SCALED  LOAD-INCREMENTS  IN 
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■6 


Figure  1.  The  residual  tangential  stress  distribution  in  an 

overstrained  tube  (b/a  =  2.0,  v  =  0.3,  u  =  0.1,  n  =  100). 
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•  2 
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Figure  2.  The  residual  radial  and  axial  stress  distribution  in 

an  overstrained  tube  (b/a  =  2.0,  v  =  0.3,  a)  =  0.1,  n  =  100). 
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